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^ ! Abstract 

cn 

■ We review some of the few available nonadiabatic calculations of 

00 

I dipole moments. We show that those carried out in a laboratory-fixed 

On ' 
O 

I able ones in a molecule-fixed reference frame. For completeness we 



set of coordinate axes are bound to fail and discuss the more reason- 



also describe the standard Born-Oppenheimer calculations of dipole 
moments. We briefly address the experimental estimation of dipole 
moments from the Stark shift of spectral lines and argue that it does 
not provide such property but a sort of energy-weighted average of 
dipole transition moments. 
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1 Introduction 



Many properties of molecular aggregates are attributed to the distribution of 
charges in the constituent individual molecules or what is commonly called 
the molecular dipole moment. This molecular property is also responsi- 
ble of the most salient features of the absorption and emission spectra of 
molecules [1-3]. Molecular dipole moments can be obtained, for example, 
from microwave spectra [1-3], or calculated from first principles [4]. 

The purpose of this review is to outline the calculation of molecular dipole 
moments by means of quantum-mechanical approaches and compare such 
theoretical results with corresponding experimental measurements. It is not 
expected to be exhaustive because of the enormous number of worthy works 
on the subject. However, we hope to succeed in giving an idea of the diffi- 
culties encountered in such endeavour. 

In Sec. [2] we consider the nonrelativistic molecular Hamiltonian and dis- 
cuss the separation of the center of mass in detail. We derive general expres- 
sions that may be useful for most of the various nonadiabatic approaches 
that appear in the current literature. In Sec. [3] we outline the dynamical 
symmetry of the molecular Hamiltonian and some of the properties of the 
molecular stationary states. In Sec. IHwe discuss the well-known Hellmann- 
Feynman theorem for variational wavefunctions because it is relevant to the 
calculation of dipole moments. Although we are mainly interested in nona- 
diabatic calculations of dipole moments, in Sec. [5] we consider the separation 
of the electronic and nuclear motions that leads to the Born-Oppenheimer 
and adiabatic approaches to the calculation of molecular properties. For 
completeness, in Sec. IHlwe outline the interaction between a system of point 
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charges and an external electric field with the purpose of defining dipole mo- 
ment, polarizabilities and other molecular properties. In Sec. [7] we briefly 
describe the models used to derive the dipole moment from the Stark ef- 
fect in rotational spectra. They will be useful for comparison between ex- 
perimental measurements and theoretical calculations. In Sec. [H] we briefly 
review the theoretical calculation of molecular dipole moments under the 
Born-Oppenheimer approximation. In Sec. [9] we discuss some of the exist- 
ing nonadiabatic calculations of dipole moments which is the main topic of 
this review. In Sec. [10] we briefly describe the more rigorous nonadiabatic 
calculations of dipole moments in the molecule-fixed reference frame. For 
simplicity we restrict ourselves to diatomic molecules because they appear to 
be the only ones accessible to current nonadiabatic approaches. In Sec. [TTlwe 
outline the application of perturbation theory to the Schrodinger equation 
for a molecule in an external electric field with the purpose of discussing a 
more rigorous connection between the measured Stark shift and the molec- 
ular dipole moment. Finally, in Sec. [12] we summarize the main conclusions 
of this review. 
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2 The molecular Hamiltonian 



In this section we consider the nonrelativistic molecular Hamiltonian as a 
system of charged point particles with only Coulomb interactions 



In this expression rrii is the mass of particle i, gj = — e or = ZjC are 
the charges of either an electron or nucleus, respectively, and Vij = |rj — rj| 
is the distance between particles i and j located at the points rj and rj, 
respectively, from the origin of the laboratory coordinate system. In the 
coordinate representation = —ihVi. 

Since the Coulomb potential V is invariant under translations f/(a)rjf/ (a)^ = 
Ti+a {U{a)VU (a)^ = V), then the eigenf unctions of the translation-invariant 
Hamiltonian operator ([1]) are not square integrable. For that reason we first 
separate the motion of the center of mass by means of a linear coordinate 
transformation 



H 



T + V, 

^ ^2 




V 




(1) 




(2) 



that leads to 




(3) 



and 



E -V? = E - + 2 E E ^v;v; + E E E ^v;v; (4) 

rrii ^-^ rrii ^-^ ^-^ rrii ■' ^-^ ^-^ ^-^ rrii ■' 

i i i i>l i j>l k>l 
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It is our purpose to keep the transformation ([2]) as general as possible so 
that it applies to all the nonadiabatic approaches discussed in this review. 
In order to uncouple one of the new coordinates r'^^ from the remaining ones 
we require that the coefficients of the transformation ([2]) satisfy 

J2^-^ = 0,j>l (5) 
^-^ rrii 

i 

The new coordinates transform under translations as follows: 

?7(a)r;.f/(a)t = r;. + a5^t,, (6) 

i 

If we require that r'^ transforms exactly as the original variables and that the 
remaining r'j are translation invariant we have 

i 

If we choose tu = ^rrii, where C, is an arbitrary real number, then Eq. ([5]) 
becomes Eq. ([7!) for j > 1. If we then substitute tu = ^rrii into Eq. with 
j = 1 we conclude that 

tii = J^,M = J2rni (8) 

i 

and r'^ results to be the well known coordinate of the center of mass of the 
molecule [5]. The choice of the coefficients of the transformation ([2]) for the 
remaining variables r'j j > 1 is arbitrary as long as they satisfy Eq. ([7j). 
Finally, the total Hamiltonian operator reads 



j>l k>l \ i * / ^ i=l j=i+l 
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where Hm is the internal or molecular Hamiltonian operator. The explicit 
form of the interparticle distances Tij in terms of the new coordinates 
may be rather cumbersome in the general case. We will consider it in the 
particular applications discussed in subsequent sections. 

For future reference it is convenient to define the center of mass and 
relative kinetic energy operators 

fcM = -I^Vf (10) 

j>l k>l \ i ' / 

respectively, so that T = Tcm + T^eu Hm = T^ei + V , and H = Tcm + Hm- 
The inverse transformation exists and gives us the old coordinates in 
terms of the new ones: 

r' = E(^"').r; (12) 

According to equations and ([7]) we have U{a.)riU{ay = (t^^)^ a+rj from 
which we conclude that 

(t-i)^^ = l, z = l,2,...,N (13) 

In order to understand the meaning of this result notice that the momentum 
conjugate to r- is given by the transformation 

p: = E(^"%p.- (14) 

so that the linear momentum of the center of mass 



p'i = Ep^- (15) 



is precisely the total linear momentum of the molecule. We also appreciate 
that TcM = p'iV(2M) and that he inverse transformation of the momenta is 
(see Eq. ©) 

i 

Eq. (IT^ is also relevant to the behaviour of the dipole moment in internal 
coordinates. We have 

f^ = '^(ltri (17) 

i 

in the laboratory-fixed frame, and 

[/(a)/if/(a)t = /i+ga, g=^gi (18) 

i 

clearly shows that /i is invariant under a translation of the origin for a neutral 
molecule q = 0. In internal coordinates we have 

/^=E^>;' ^; = E^^(t"')., (19) 

j i 

and f/(a)/if/(a)'l" = /i+g^a that is consistent with Eq. f|T8l) because q'l = qhy 
virtue of Eq. f[T^ . 

We have kept the transformation ([2]) as general as possible in order to 
have a suitable expression of the molecular Hamiltonian that applies to all 
the nonadiabatic approaches discussed in this review. In what follows we 
illustrate some particular ways of determining the remaining transformation 
coefficients. More precisely, we will choose a point in the molecule as the 
origin of the new coordinate axis and refer the positions of — 1 particles 
to it (remember that three coordinates have been reserved for the location 
of the center of mass). 
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As an example, consider an arbitrary set of particle labels / and choose 
' ' Mj 



r', = r.- - -j^ ^ niiTi, j > 1, = ^ (20) 



~ Mr 

iel i&I 

That is to say: 



rrii 



1 if i G / 



tji = 6ji - -Ti-Sii, j > 1, Sii = < (21) 
^^■f I otherwise 

These coefficients tji already satisfy Eq. ([7]). If / is the set of nuclear labels 
then we refer the positions r^, j > 1 to the nuclear center of mass. If the 
set I contains only one nuclear label we refer the positions of the remaining 
particles to that particular nucleus. Different authors have already chosen 
one or the other coordinate origin as shown in sections [9] and [TOl 
If we take into account that 

j^t^^S^_l ^ ^ 1 

rrii rrii Mj Mj 

i ■' 

then we realize that the transformation fl20|) uncouples the coordinates of the 
particles that belong to / from the remaining ones in the kinetic operator: 

' ^-^ - ^ a j,k e I 

if j G / and k ^ I oi j ^ I and /c G / (23) 



Mi 

The equations derived in this section are sufficiently general to cover all 
the strategies commonly followed in the separation of the center of mass prior 
to the solution of the Schrodinger equation for atoms and molecules. Other 
general expressions for the translation-free internal coordinates have already 
been discussed by Sutcliffe and Woolley [6-9]. It is a most important topic 
in the nonadiabatic quantum-mechanical approach to atoms and molecules. 
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3 Symmetry 



In this section we discuss the symmetry of a molecular system that is deter- 
mined by all the operators that commute with its Hamiltonian. For example, 
the total Hamiltonian operator commutes with the total linear momentum 

N 

[H,p] = 0, p = 5^p, =p; (24) 
and also with the total angular momentum 

TV 

[H,J]=0,3 = Y.h (25) 

i=i 

It follows from equations (fT2l) and (fT6!) that 

J = ^f, X p, = J^f^ X p: = f; X p; + J' (26) 

i i 

where J' is the internal or molecular angular momentum. It is clear that 

N 

[^M,J']=0, 3' = Y,i[^t>[ (27) 

The Hamiltonian operator is invariant under permutation of identical 
particles. We may formally write this symmetry as 

VHV-^ = H ^[H,V]=0 (28) 

where V stands for any permutation of electrons or identical nuclei. Since 
TcM is invariant under permutation of identical particles we conclude that 

VHmV = Hm [Hm, V] = (29) 
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Besides, the molecular Hamiltonian operator is also invariant under in- 
version 

^Hm^~^ ^ [Hm, ^] = (30) 

where the inversion operator $J produces the transformations ^rj^~^ — —rj 
and ^pj^-^ = -pj. Also notice that [J',^] = 0, [J',V] = and [Q,V] = 
so that we can choose a stationary state ^ to be simultaneous eigenfunction 
of such set of commuting operators: 

j'H = /iV(J' + 1)^, J' = 0, 1, . . . 

j;* = hM'j^, m; = o,±i,...,±j' 

= ±^ (31) 

where u = 1 or cr = —1 for identical bosons or fermions, respectively. A some- 
what more rigorous discussion of the kinematics and dynamics of molecules 
was provided, for example, by WoUey [10]. The simplified discussion in this 
section is sufficient for present purposes. 

4 The Hellmann— Feynman theorem 

In this section we derive the Hellmann-Feynman theorem for optimized varia- 
tional wavefunctions because it is most important for the accurate calculation 
of dipole moments [4, 11]. 

Consider a trial function $ and the variational energy W given by 

W {^\ ^) = {^\ H \^) (32) 
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An arbitrary variation 5$ leads to 

5W ($1 ^) = {5^H -W \^) + {5^H -W 1$)* (33) 

were the asterisk denotes complex conjugation. The right-hand side of this 
equation vanishes when 5W = 0; that is to say, when W is stationary with 
respect to the infinitesimal change 5$ in the trial function. The variational 
or optimal trial function satisfies this condition. 

If the Hamiltonian operator depends on a parameter A, which may be, for 
example, a mass, charge, force constant, the strength of a external field, etc., 
then the variational function will also depend on it. Therefore, differentiation 
of Eq. (1521) with respect to A leads to 
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H-W\^Y (34) 



It is clear from this equation that if the set of variations 5$ includes {d^/dX) 5\ 
then the optimized trial function $ will satisfy the well known Hellmann- 
Feynman theorem [12,13]: 



where (A) = ($|A|<I>) /($|$). 




' • (35) 



5 The Born— Oppenheimer Approximation 

Since the vast majority of quantum-mechanical studies of molecular proper- 
ties are based on the Born-Oppenheimer (BO) approximation [14] we discuss 
it here with some detail, in spite of the fact that we are more interested in 
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non BO approaches. Before proceeding we mention the curious fact that the 
modern and most useful form of that approach appears to have been first 
developed by Slater [15] in his study of the Helium atom as early as 1927. 
However, we will adopt the common practice to attribute it to Born and Op- 
penheimer. In this section we begin with the traditional procedure proposed 
by Slater [15] and Born and Huang [14] and then mention other authors' 
approaches and criticisms. 

For simplicity and clarity we write the Schrodinger equation for a molecule 

as 

^^,(r^r") = E,^„(r^r'^) (36) 

where and r*^ denote all the electronic and nuclear coordinates, respec- 
tively, and the subscript a stands for the collection of all the quantum num- 
bers necessary for the description of the molecular state. We use the symbols 
Te, T„, V"ee,V"„e, Km to indicate the kinetic energy operators for the electrons 
and nuclei, and the electrons-electrons, nuclei-electrons, and nuclei-nuclei 
Coulomb interactions, respectively: 

H = fn + He + Kn, He = t + Ke + Ke (37) 

The BO approximation is based on the ansatz 

vI/,(r^r'^) = ^X«,(r")^,(r^r") (38) 

j 

and we assume that {ipj \ fk)e = ^jk, where the subscript e indicates integra- 
tion over electronic coordinates. It follows from equations fl36l) . fl37j) and fl38l) 
that 

{(fk\ H |^c.)e = TnXak + (Ukj + Wkj) Xaj = E^Xak (39) 

j 
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where 

Ukj = {ipk\ He \(pj)^ + VnJkj, Wkj = {ifikl [fn, ^j] ) (40) 

/ e 

We have introduced the formal commutator [T„, ipj] with the sole purpose of 
making the resulting equation more compact. 

In the standard BO approximation we choose the functions (pj to be 
eigenfunctions of the electronic Hamiltonian Hg 

4¥',(r^r") = E,,(r")¥,,(r^r") (41) 

so that 

f/fc,(r") = f/,(r")5,,, f/,(r") = E,,(r") + Kn(r'^) (42) 

In the first approximation we assume the couplings Wkj between "electronic 
states" to be small and keep just one term for the ground state 

vP,(r^r")^X.o(r")¥^o(r^O (43) 

and 

f„ + f/o(r")] xaoin = E^'^xaoin (44) 

The addition of the diagonal term Wqq to the nuclear operator in this equation 
gives rise to the adiabatic approximation. 

Sutcliffe and Wolley have criticized the main assumptions of this approach 
in several papers [7-10,16]. For example, they pointed out that the functions 
in the expansion in the right-hand side of Eq. (1551) are assumed to be square 
integrable, whereas the left-hand side is not because we have not removed 
the motion of the center of mass. Such criticism, which also applies to the 
approximate and widely used ansatz fH3|) . does not appear to be justified 
because the nuclear functions Xaj are not square integrable. According to 
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those authors another drawback of the BO approximation is that it does 
not take into account the permutational symmetry of identical nuclei that 
are treated as distinguishable particles clamped in space in order to define 
a framework geometry unambiguously. However, in principle we can intro- 
duce such permutational symmetry in the approximate wavefunction fH51) 
by means of appropriate projection operators like those mentioned below in 
Sec. m In our opinion the Born-Oppenheimer approximation is a consistent 
way of deriving a suitable approximate solution to the Schrodinger equation. 

It is interesting how some popular good books on quantum chemistry 
introduce this subject. For example, Szabo and Ostlung [4] state that "Our 
discussion of this approximation is qualitative. The quantitative aspects 
of this approximation are clearly discussed by Sutcliffe" and they give a 
reference and simply write the single product (H3|) . On the other hand. 
Pilar [17], referring to a Hamiltonian operator like ([1]), says: "where it is 
assumed that all nuclear and electronic coordinates have been referred to 
the center of mass of the system". However, that Hamiltonian operator 
exhibits all the electronic and nuclear coordinates and no coupling terms 
coming from the separation of the center of mass. Therefore, it cannot be 
the internal Hamiltonian operator as we have already seen above. These are 
just two examples of how most authors blindly accept the clamped nuclei 
approximation without further analysis. 

In some cases, mostly for diatomic molecules, the separation of the motion 
of the center of mass has been carried out rigorously prior to the application 
of the BO approximation [12]. Sutcliffe and Woolley have already discussed 
this issue in detail in several papers [7-10,16]. 
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Although it is not customary to remove the motion of the center of mass 
of the whole molecule before the application of the BO approximation, it 
appears to be common practice to separate the nuclear center of mass at the 
second stage (jS]). Since the nuclei are much heavier than the electrons we 
may assume that the error is small. A rigorous discussion of this issue as 
well as a simple exactly solvable example are given elsewhere [18,19]. 

A most important by-product of the clamped-nuclei approximation is 
that it enables one to introduce the familiar (classical) chemical concept of 
molecular structure into the quantum-mechanical calculations. The molec- 
ular geometry is given by the equilibrium nuclear configuration r^^ at the 
minimum of [/(r"): 



At this point it is worth noticing that U (r") would not appear in a straight- 
forward rigorous solution of the Schrodinger equation fl36|) . Sutcliffe and 
Woolley [6-10, 16, 20] have also pointed out to the difficulty of discussing 
molecular structure without the BO approximation (or any of its variants). 

6 Interaction between a molecule and an ex- 
ternal field 

In order to make the present review sufficiently self-contained we briefly 
develop the main equations for the interaction between a molecule and an 
external potential $(r). The energy of that interaction is given by [21] 



V^(r'^)lr- 



=r; 







(45) 




(46) 
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where p(r) is the molecular charge density. If the potential varies slowly in 
the region where the charge density is nonzero we can expand it in a Taylor 
series about the coordinate origin located somewhere in the molecule: 

*(r) = *(0) + J: u^(0) + ^ ^ ^ ™|±(o) + . . . (47) 

U U V 

where u,v — x,y,z. Since V^$(r) = we subtract r^V^$(0)/6 from this 

equation and rewrite the result as 

rl<b 1 Ff'^ 

$(r) = $(0) + E - a^(0) + 6 E E (3- - r'^u.) ^^(0) + • • • (48) 

U U V 

If we now take into account that the external field is given by F(r) = — V$(r) 
we have 

1 r) F 

$(r) = m - r • F(0) - ^ E E (3^^ " ^'^™) ^(0) + " " " (^9) 

U V 

For a set of point charges Qi located at r^, i = 1, 2, . . . , we have 

N 

p(r) = E ^^'^(^ - (50) 

i=l 

and the interaction energy becomes 

1 rlW 

w = qm - • F(o) - g E E ^-^(0) + ■ ■ ■ (51) 

U V 

where the net charge q, the dipole moment /i, and the quadrupole moment 
Q are given by 

N 

q = E^^ 



i=l 

N 

i=l 
N 

Quv = E (3^»^i ^ 'T'l^uv) (52) 



i=l 
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If the applied field is uniform in the region of interest, then all the terms 
beyond the first two ones in the right-hand side of Eq. flSTj) vanish. 

In one of his well known discussions of long-range intermolecular forces 
Buckingham [22] considers a Hamiltonian operator of the form 

U V 

where = Q/2 and F^v = {dFu/dv){0). Thus, the perturbation expansion 
for the energy of the molecule in the external field results to be [22] 

U U V U V 

(54) 

where 

/.f = (v^(°)|/i.|vl/W> 

el°i = (Vl/W|0™|vl/W> (55) 

are the permanent dipole and quadrupole moments, respectively. Bucking- 
ham [22] does not write the Hamiltonian operator H^^^ explicitly, but he 
refers to the energy of "separate molecules for fixed molecular positions and 
orientations". Therefore one may assume that he probably means the BO 
Hamiltonian operator. On the other hand. Bishop [23] considers that Eq. flMl) 
can be used irrespective of whether one is considering the electronic or the 
total molecular energy. He explicitly indicates that qi is an element of charge 
at the point rj relative to an origin fixed at some point in the molecule. 
As shown below in sections M and [101 this choice of reference frame is not 
convenient for the nonadiabatic calculation of the dipole moment. 
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In a most interesting and comprehensive review of the electric moments 
of molecules Buckingham [24] argues, by means of simple and rigorous sym- 
metry arguments, that a diatomic molecule in some stationary states does 
not possess a dipole moment. We will come back to this point in sections [9] 
andM 

7 Experimental Measurements of Dipole Mo- 
ments 

Typical experimental measurements of dipole moments are based on mi- 
crowave spectra, and almost invariably on the Stark effect and the model of 
a rigid (or almost rigid) rotating dipole. If Hrot is the hamiltonian of a rigid 
rotator and — /i ■ F is the interaction between the molecular dipole moment 
and the uniform electric field (as shown above in Sec. EJ), then the Stark 
rotational energies are given by 

Hrot{FW"' = ^^°*(F)^™*, Hrot{F) = HrotiO) - /iFcos(^) (56) 

where 6 is the angle between the external field F and the dipole moment /i. 

The external electric field produces both shift and splitting of the molec- 
ular rotational energies. The magnitude of the shift of the spectral lines 
z/(F) = AE^°^{F) /{he) changes with the field. Rayleigh-Schrodinger pertur- 
bation theory provides analytical expressions for the Stark shifts in the form 
of a power series of the field intensity: 

Au{F) = u{F) - i>(0) = u^^^^iF + i>^^\fiFf + ... (57) 
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where the coefficients u^^^ are known functions of the rotational quantum 
numbers [1-3]. If we measure the Stark shifts for known values of the field 
intensity and then fit selected experimental data to a polynomial function of 
the field we can obtain the dipole moment from the polynomial coefficients. 

In order to account for the hyperfine structure of the spectra one should 
add the terms arising from the spin-rotational interaction, the quadrupole 
interaction of the nuclei, and the spin-spin magnetic interactions [1-3]. 

Molecular beam electric resonance experiments are also based on the same 
model of rotating dipole [3] . 

There have been many experimental studies of a wide variety of molecules. 
We restrict to the simplest ones that are accessible to existing nonadiabatic 
approaches, such as, obviously, diatomics [25-28]. In this case u'^^^ = and 
the Stark effect is quadratic. It is worth paying attention to the discrepancy 
in notation and presentation of the empirical models. For example, the model 
Hamiltonian may be either an actual operator [27] or a scalar function of 
the quantum numbers [28]. The interaction between the dipole and the 
field may also be written in somewhat different ways [27, 28] or the effective 
Hamiltonian may omit the rotational kinetic energy [26]. 

Although symmetric tops appear to be beyond present nonadiabatic treat- 
ments, we quote them here as another example of the use of the model of a 
rotating dipole outhned above [29,30]. In this case i>*^^) = when K — and 
u^^^ 7^ when K ^ 0, where K is the quantum number for the projection of 
the angular momentum along the symmetry axis. 

It is clear from the discussion above that the accuracy and reliability of 
the experimental determination of the molecular dipole moment does not 
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depend only on the accuracy of the measured Stark hne shifts but also on 
the validity of the theoretical, semi-empirical model of the molecule as a 
rotating electric dipole. We will come back on this important point later in 
Sec. EH 



8 Born— Oppenheimer Calculations of Dipole 
Moments 

Although we are mainly interested in the nonadiabatic calculation of dipole 
moments it is worth comparing it with the more popular BO approach to 
the problem. For this reason, in this section we outline the latter. If we are 
able to solve Eq. fHT!) for an appropriate set of nuclear configurations and 
determine the equilibrium geometry of the molecule given by Eq. (H5l) . we 
can then calculate the dipole moment as follows: 



for a molecule with N^, electrons and A^^ nuclei [4]. In this equation dr^ 
denotes the volume element for all the electronic coordinates. 

An alternative approach consists in solving the electronic BO equation 
for the molecule in an electric field 

ife(F) =i^e(0)+eF-^r^^ (59) 

i=i 

and then calculate the electronic part of the dipole moment as 
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where Ee{r^g,F) is the lowest eigenvalue of the electronic operator fl59l) for 
the equilibrium geometry. 

Szabo and Ostlund [4] discuss the reasons for the noticeable discrepancy 
between both approaches. It is known to be due to the fact that approximate 
wavefunctions which are not well optimized fail to satisfy the Hellmann- 
Feynman theorem [11] already discussed in Sec. |H Since the formally correct 
definition of properties like the dipole moment is as a response function to 
an external field (see Sec. [H]) Swanton et al [11] proposed a calculation based 
on the following expression 



i^e(0)-i?e(r:,,0)|$(0)) (61) 



in such cases where the variational function is not fully optimized and one 
does not have appropriate analytical expressions for the derivatives {dE^j dFy) f=o- 
This equation is a particular case of (IMI) when the electronic approximate 
wavefunction is real and normalized to unity ($(0)| $(0)) = 1. 

The nuclear configuration given by the set of equilibrium coordinates r^^ 
determines what we usually call the geometry of the molecule and thereby 
the orientation of the dipole moment. For example, we know that the dipole 
moment is directed along the molecular axis in a linear molecule without 
an inversion center or along the symmetry axis in a symmetric top. The 
analysis is not so simple in the case of the nonadiabatic calculations that we 
will discuss in Sec. M 

According to Cade and Huo [31] the dipole moment calculated by means 
of this quantum-mechanical approach for just one internuclear distance (say, 
the theoretical or experimental R^) is not strictly comparable to the exper- 
imental one that commonly corresponds to a particular vibrational state or 
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an average over a set of vibrational states. In spite of this apparent defi- 
ciency of the BO approach, it is worth noticing that when Wharton et al [25] 
determined the dipole moment of LiH experimentally to be = 5.9D there 
were as many as eight previous reliable quantum-mechanical calculations 
that agreed with it to ±0.3D [32]. 

9 Nonadiabatic Calculations of Dipole Mo- 
ments 

Most theoretical calculations of dipole moments are based on the BO approx- 
imation for several reasons: first, and most importantly, non BO calculations 
require far more computer time, second, most theoretical chemists are un- 
willing to go beyond the BO approximation that they assume to know well, 
third, non BO calculations give rise to some additional theoretical difficulties 
as we will show below. It is therefore not surprising that non BO calcula- 
tions of dipole moments are scarce and few. We describe some of then in this 
section and in Sec. [TOl 

Tachikawa and Osamura [33] proposed a dynamic extended molecular 
orbital (DEMO) approach based on SCF wavefunctions of the form 

^SCF^Yl^I (g2) 

I 

where each $^ is a function of the coordinates of a set of identical particles 
with the appropriate permutation symmetry. These functions are expressed 
in terms of generalized molecular orbitals 0j that are linear combinations of 
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floating Gaussians xi- 

^] = E "^nXl (63) 

r 

It is worth noticing that Tachikawa and Osamura[TO00] did not separate 
the motion of the center of mass, and this omission will give rise to consider- 
able errors if the SCF orbitals depend on laboratory-fixed coordinates [34]. 
Another limitation of this approach is that the SCF wavefunction does not 
take into account particle correlation that may be quite strong between nu- 
clei [35,36]. 

If the variational approach were based on a trial function (/^(rg, rg, . . . , r'^) 
of just internal, translation-free, coordinates r^- (see Sec. [2]), it would not be 
necessary to separate the motion of the center of mass explicitly because 
{ip\H\ip) = {ip\ Hm \ f) ■ But it is not the case of the SCF wavefunction 
m so that we have E^^^ = (^^^^| H |^^^^) = (^^^^| Tcm l^'''''") + 
(^^scF^ fj^^ ^x^scF-^ > ^vl/^c'^l Hm I'if^'^^). Therefore, the estimated energy 
will always be worse than when the trial function depends only on internal 
coordinates, even though the SCF wavefunction may satisfy the virial theo- 
rem 2 (^^^^1 f 1^-5^^^) = - (^^JSCF^ Y [13,33]. The reader may find 
a more detailed discussion of this issue elsewhere [34]. 

Tachikawa and Osamura [33] calculated the dipole moments of the ™H"H 
and ™Li"H isotopomer series, but, unfortunately, they did not show the ex- 
pression that they used. This issue is not a minor one as discussed in what 
follows. 

Before proceeding further, it is convenient to discuss the failure of the 
naive approach to the nonadiabatic calculation of dipole moments. In quan- 
tum mechanics one obtains the average of an observable O as the expectation 
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value O of the corresponding operator O. In Sec. Owe showed that 
53\E' = ±\& because the molecular Hamiltonian is invariant under inversion. 
Since 53/153^^ = —jl we conclude that 

^ = {m\jl\m)={] (64) 

for any nondegenerate molecular state ^ . This result that is known since 
long ago [24] applies to any molecule in its ground state and renders fruitless 
the calculation of its dipole moment as a straightforward expectation value in 
a set of axes parallel to the laboratory one (like the one discussed in Sec. |2]). 

One of the problems that arises from the separation of the center of mass 
discussed in Sec. H] is that the transformation fll2p can make the Coulomb 
potential rather messy. One way of keeping a simple form of the potential- 
energy function is to choose one of the particles as coordinate origin. From 
a practical point of view it appears to be convenient to choose the heaviest 
nucleus for that purpose [35]. Thus, the transformation 

N 
i=l 

r;. = r,-ri, j = 2,3,...,iV (65) 

leads to 

i j>l k>l i=l j=i+l ■' 

where mi is the mass of the heaviest nucleus and ri its location in the 
laboratory reference frame. The transformed Coulomb potential does not 
make the calculation of matrix elements unnecessarily complicated as fol- 
lows from the fact that |rj — ri | = |r^| , i = 2, 3, . . . , and |r,j — | = — r^- 1 , 
^,j = 2,3,...,iV. 
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As an example we consider a four-particle molecule. The transformation 
(|65|1 is given by the matrix 



(67) 





mi 
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m4 
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-1 
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/ 



with inverse 



.-1 
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m2 



1714 



M 
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m2 
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m4 


M 
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m2 


1 rn:i 


m4 
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^ M 
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m2 


ms 


1 m4 


M 


M 





(68) 



This transformation applies, for example, to the H2 isotopomer series [37]. 
The only members that exhibit dipole moments are HD, HT and DT. In this 
case we choose the labels 1,2,3,4 to denote the heaviest nucleus, the lightest 
one, and the two electrons, respectively. 

The fact that nii > 1112 gives rise to a charge asymmetry and a small 
dipole moment. Obviously, the straightforward BO approximation cannot 
account for it because the electronic Hamiltonian He does not depend on 
the nuclear masses (see Sec. E]) and the resulting electronic charge density is 
exactly the same for all isotopomers. Therefore, in order to obtain the dipole 
moment of diatomic molecules of the form '"A'^A, one has to take into account 
nonadiabatic corrections [38-40]. An alternative approach is based on the 
fact that there is no unique way of implementing the BO approximation. In 
fact, appropriate canonical transformations of the coordinates prior to the 
application of the BO approximation may force the required asymmetry and 
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provide a suitable way of calculating the dipole moment of "^A"A molecules 
[41,42]. 

Cafiero and Adamowicz [35, 37] calculated nonadiabatic dipole moments 
for some small diatomic molecules; in what follows we outline the variational 
method proposed by those authors for a diatomic molecule with N — 2 elec- 
trons. The core of the approach is a basis set of floating s-type explicitly 
correlated Gaussian functions of the form 

^fe(r) = exp [-(X - Sfc^) • Afe • (X - Sfc^)* - (Y - s^y) • • (Y - Sfe^)* 
-(Z-Sfe,)-A,-(Z-Sfc,)*] (69) 

where Ak is a N' x N' {N' = N — 1) symmetric matrix, X, Y, and Z are 
1 X A^' matrices of the form X = (xg, x'^, . . . , x'j^), is a 1 x A^' matrix that 
determines the location of the center of the Gaussian in space and t stands 
for transpose. In order to assure that the Gaussians are square integrable 
they chose A^ to be of Cholesky factored form A^ = • L^., where is a 
N' X N' lower triangular matrix. 

In order to have the correct permutation symmetry of identical particles 
they resorted to appropriate projection operators of the form [35] 

E = l[Ei (70) 

i 

and constructed the variational function 

m 

^ = J2ckEgk{r) (71) 

k=l 

SO that the permutation of any pair of identical particles leads to either \1' or 
— ^ if they are bosons or fermions, respectively. Then, they minimized the 
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variational energy 

E = min V' '/ 72 

Notice that the matrix has N'{N' + l)/2 independent variational param- 
eters and each contributes with A^'; therefore, there are N'{N' + l)/2 + 
3A^' + 1 adjustable variational parameters for every basis function. In a most 
comprehensive review Bubin et al [43] discussed this type of calculation in 
detail. 

Since the straightforward expectation value of the dipole-moment oper- 
ator /t will not produce any physically meaningful result, as discussed above 
in this section, Cafiero and Adamowicz [35-37] resorted to an alternative 
approach based on the energy of the molecule in an external electric field F, 
given by the Hamiltonian operator: 

H{F) = HM-F-fi (73) 

They fitted energy values to a polynomial function of the field [35-37]: 

E{F,) = E{0) - - - . . . (74) 

and obtained the dipole moment from the linear term. The dipole moments 
of HD, HT, LiH and LiD calculated in this way proved to be very accurate 
and in remarkable agreement with available experimental values [35,37]. In 
particular, the rate of convergence of the theoretical results for the LiH [35] 
towards the corresponding experimental dipole moment [26] is astonishing. 
Table [1] shows some dipole moments calculated by Cafiero and Adamowicz 
[35,37] and other authors [33,44,45] as well as the corresponding experimental 
values [26,46]. 
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However, if the approximate variational function (171 p were fully opti- 
mized, then it should satisfy the Hellmann-Feynman theorem discussed in 
Sec. m and in that case the only possible result would be 

dE{F, 



= (^|/i,|^) = (75) 

The only way that Cafiero and Adamowicz [35, 37] could obtain a nonzero 
dipole moment by means of this approach is that their variational wavefunc- 
tion did not become spherically symmetric when — > 0. In other words, 
that their variational ansatz was not sufficiently accurate at small values of 
the field strength where it would reveal the permanent molecular dipole mo- 
ment. In that case, it is not clear how those authors [35,37] obtained such 
remarkable agreement between their theoretical results and the experimental 
values of the dipole moment [26]. A more detailed discussion of this baf- 
fling agreement has been published elsewhere [47, 48] . A most speculative 
explanation is that a biased placement of the centers of the floating Gaus- 
sians somehow mimicked the use of a body-flxed set of axes (see Sec. [TU]) . 
However, the authors never mentioned this possibility. Curiously, in a later 
paper the authors state that "This spherical symmetry for the ground-state 
wavefunction implies several things about molecules that may go against 
common chemical intuition. First of all, no molecule in the ground state will 
have a dipole moment, just as atoms do not. Similarly, the molecule will 
have only one unique polarizability, an isotropic polarizability. The current 
authors have presented several papers which discuss these phenomena" [49]. 
However, they did not explain how they had obtained the dipole moments 
in two of those earlier papers [35,37]. Besides, in Sec. [10] we will show that 
the statement "no molecule in the ground state will have a dipole moment, 
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just as atoms do not" is false. 

It would be interesting to investigate to which extent the variational 
ansatz flTTl) optimized by Cafiero and Adamowicz [35,37] satisfies the Hellmann- 
Feynman theorem that in the present case gives us an exact theoretical rela- 
tion between the molecule's response to the field and the dipole moment. Ex- 
act relationships like those given by the hypervirial and Hellmann-Feynman 
theorems are useful to determine the accuracy of approximate wavefunc- 
tions [12,13]. The questionable success of the method devised by those au- 
thors is obviously based on the fioating nature of the Gaussian functions 
that one can place conveniently to get the desired result [35,37]. Notice that 
Bubin et al [50] resorted to one-center Gaussian functions of the form 

0,(r) = exp [-X • A, • X* - Y • ■ Y* - Z ■ ■ Z*] (76) 

to determine the charge asymmetry of the rotationless states of the HD 
molecule. Clearly, this kind of basis functions cannot be placed at will to force 
an axial symmetry and a nonzero dipole moment. It is for this reason that 
we mentioned above in this section that it is unfortunate that Tachikawa and 
Osamura [33] did not show the expression that they used for the calculation 
of the dipole moment of the ™H"H and '"Li'^H molecules. 

The nonadiabatic calculations of molecular dipole moments described in 
this section are bound to fail because they are based on a set of axes parallel to 
the laboratory one. Therefore, the expectation value of the dipole-moment 
operator vanishes for any molecule in a nondegenerate state. In Sec. [10] 
we discuss more judicious calculations based on molecule-fixed coordinate 
systems. 
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10 Molecule— Fixed Coordinate System 



It is convenient to discuss the motion of a system of particles in space by 
means of three sets of axes. One set of axis remains fixed somewhere in 
the laboratory. A second set of axes with its origin fixed on the system 
center of mass and parallel to the laboratory one. It was discussed in Sec. [2] 
with the purpose of separating the motion of the entire system as a point 
particle with mass equal to the total mass of the system. A third set of 
axes with origin on the center of mass and somehow completely fastened to 
the system enables us to describe its rotational motion. It is not difficult to 
define such set of axes for a rigid body, but it is not so obvious when the 
distance between the system particles change rather arbitrarily. It is clear 
that we should define any intrinsic molecular property, like the electric dipole 
moment, moment of inertia, etc, with respect to this body-fixed reference 
frame. If we use the second set of axes we expect that the expectation value 
of the dipole-moment operator for the molecule in a nondegenerate state 
vanishes as discussed in Sec. M This result is an obvious consequence of the 
average over the angular degrees of freedom of the whole system. The proper 
use of the body-fixed set of axes was clearly addressed by Blinder [38, 39] 
in his studies on the HD molecule, and Sutcliffe [6-8, 51] discussed the more 
general case of polyatomic molecules. For simplicity, we restrict ourselves to 
Blinder's proposal for diatomic molecules [38, 39] in what follows. 

For generality we first consider a diatomic molecule with the nuclei located 
at ri and r2 and the electrons at r2, rs, . . . , r^v. Following Blinder [38] we use 
relative coordinates for the nuclei and refer the electron coordinates to the 
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midpoint between the nuclei according to 

i 

r's = ra - ri 

< = r.-^(ri + r2), i>2 (77) 
that is a particular case of Eq. ([2]) with 

TTLj 

hj = Sj2 - Sji 

U, = 5., -i (5, 1 + 5,2), J>2 (78) 

The resulting molecular Hamiltonian operator in the coordinate representa- 
tion will be 

^"-^^-^5:v?-^ Ev; +^v'.Ev;.v^ (79) 

" " j=3 * \i=3 / " i=3 

where me is the electronic mass, rris = mima/ (mi+ma) and = mim2/(mi — 
ma). These coordinates are most convenient to calculate the dipole moment 
of diatomic molecules of the form '"A'^A because the fourth term in the right- 
hand side of the internal Hamiltonian fl79|) vanishes when mi = ma and can 
therefore be treated as a perturbation that converts the symmetric case to 
an asymmetric one. 

In order to place the molecule-fixed set of axes we take into account the 
unit vectors generated by 

R = ra — Ti = i?(cos sin 96^ + sin sin 6 ey + cos 6 e^) (80) 
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where e-^, ey and are the space-fixed orthonormal Cartesian vectors (sec- 
ond set of axes). Notice that we renamed the vector r'g in order to match 

Bhnder's notation [38]. 

We next define the body-fixed orthonormal vectors 



eg 
eel 



I 



OR 



and the new particle coordinates with respect to them 



(81) 



(82) 



The transformation between these two sets of coordinates 



c 



is obviously given by the orthogonal matrix 



^ e • e' 



c 



^X ' ' ^y ' ^y 



\ 



( 



I 



COS COS Q sin cos Q 



— sin cos 
y cos sin ^ sin sin ^ 



(83) 



-sin^ 


cos^ 



(84) 



We do not show the explicit form of the Hamiltonian operator in this molecule- 
fixed reference frame because it is not necessary for the present discussion. 
Blinder [38] derived it for the HD molecule. 
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The new nuclear coordinates will be 



Ri = ri-- 



R2 = r2 - - 




) 



(85) 



and the electronic ones are simply given by r^', i > 2. Assuming that both 
nuclei have identical charges then the classical dipole moment in the body- 
fixed set of axes is given by a purely electronic contribution 



For the particular case of HD {N — 4:) the expressions above agree with 
those developed by Blinder [38] , except for the different notation and labeUing 
of nuclear and electronic coordinates. They allow one to calculate the dipole 
moment as the expectation value of the corresponding operator by means of 
an eigenstate of the Hamiltonian operator in the molecule-fixed set of axes. 

The effect of the inversion operator on the rotation angles is given by 
9 ^ TT — 9 and (p ^ (f) + n. Therefore, the first row of the matrix C remains 
unchanged and the other two ones change sign. For that reason the inversion 
operation (a;^,y^,^Q —>■ {—x'^, —y[, —z[) in the laboratory-fixed frame results 
in {x'l,y'l,z'l) — > {—x",y",z") in the body-fixed one, and, consequently, we 
do not expect that (ji") vanishes because of inversion symmetry. However, 
in the case of identical nuclei rrii = m2 we cannot have a net dipole moment 
along the intcrmiclcar axis because the additional permutational symmetry 
leads to {fi'^) — 0. In order to explain the occurrence of a dipole moment in 
a diatomic molecule of the form "*A"A we resort to perturbation theory and 



N 




(86) 



i=3 
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write 



M — 



Hq + \H' , where 




2me 



8m 



'S 



( v'3 + v;)' + V 



(87) 




from the fact that A = when mi = m2- 

We can thus expand the eigenfunction in a A-power series \1/ = + 
\[^(^^ . . . that leads to a similar expansion for the dipole moment Hz = /ii^^A + 
. . . where 



provided that the eigenfunction is chosen to be real. This simple argument 
shows how the different nuclear masses produce a charge asymmetry and a 
net dipole moment. 

Blinder [38] estimated the dipole moment of HD by means of a rather 
more complicated perturbation approach and later Kolos and Wolniewicz 
[44,52] and Wolniewicz [53,54] improved Blinder's calculation by means of a 
variational-perturbation method based on Eq. [891 

It is clear that any rigorous nonadiabatic calculation of the molecular 
dipole moment should be carried out in the molecule-fixed set of axes. Con- 
sequently, those approaches outlined above in Sec. M are unconvincing (to 
say the least). Cafiero and Adamowicz [35,37] must have placed the floating 
Gaussians in a convenient way to obtain nonzero dipole moments in good 
agreement with the experimental ones and Tachikawa and Osamura [33] did 




(89) 
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not explain how they obtained their results. Besides, the latter authors even 
forgot to remove the motion of the center of mass. 



11 Perturbation theory for the Stark Shift 

As we outline above in Sec. [3, the experimental determination of the molecu- 
lar dipole moment relies on the validity of the model of a rotating quasi-rigid 
polar body. The procedure consists of fitting a polynomial function of the 
field strength to the observed Stark-shift lines. In principle, we should derive 
more rigorous theoretical expressions for those line shifts by means of per- 
turbation theory and the actual quantum-mechanical Hamiltonian operator 
for a molecule in an electric field. 

Instead of the rigid rotator model we should choose -f^o = Hm the Hamil- 
tonian operator for the isolated molecule and the perturbation H' = —F.p,. 
If we set F along the z axis in the laboratory frame (F = Fbz), then 
H' = —Ffiz- In this way, perturbation theory gives us the well-known 
quantum-mechanical expressions for the Stark shifts 



where n is a collection of quantum numbers that completely specifies a given 
(nonadiabatic) molecular state. For the ground state we have 



(90) 




(91) 
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ement ( \I/q 



(0) 



where we clearly appreciate that the field-reduced splitting AE/F"^ does not 
give us what we may call the square of the dipole moment but a kind of 
energy-weighted average over those molecular states with nonzero matrix el- 
ftz "^m^^, including the continuum part of the spectrum. This 
important issue was already discussed by Brieger et al [55] and Brieger [56] 
more than twenty years ago in their BO study of the Stark effect of heteronu- 
clear diatomic molecules in states. 

In the molecule-fixed reference frame we have 

= + IJ^'y^y + (92) 

and a similar expression for F. Brieger [56] argued that the spherical vector 
components 

<i = T^K±«;) (93) 

are more convenient for the calculation of the matrix elements in Eq. ( 19111 . 
It is not difficult to prove that 

/x.F = 5^(-i)>;f; (94) 

p 

However, since the laboratory-fixed and molecule-fixed are the natural frames 
for the external field and the molecule, respectively, then Brieger [56] chose 
F = Fez in the former and fi in the latter as in Eq. (!92i) . Thus, in the simple 
notation of Sec. [TOl we have 

= — sm6 fi'^ + COS n'^= ' {l^'+i ~ /^-i) + COS /Iq (95) 

v2 
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In this way, Brieger [56] showed that the perturbation H' connects the 
ground electronic state with excited and ^11 ones (in general, those 
with electronic angular-momentum quantum numbers A and Ail). Under 
the BO approximation Brieger et al [55] and Brieger [56] separated the Stark 
shift f pT]) into four contributions: (I) coupling of rotational states within the 
same vibronic one, (II) coupling of vibrorotational states within the same 

electronic state, (III) coupling of vibrorotational states between two 
electronic ones, and (IV) coupling of vibrorotational states between the given 

and ones. The first three contributions are due to /io and the fourth 
one to /i±i [55, 56]. 

As far as we know there is no nonadiabatic calculation of the Stark shift 
by means of equation fPTl) . A straightforward comparison of BO and non BO 
calculations may lead to some difficulties as mentioned by Wolniewicz [54]: 
"Since the familiar classification of electronic states of diatomic molecules is 
based on the Born-Oppenheimer approximation, some difficulties arise if one 
tries to use the standard nomenclature to describe nonadiabatic functions" 

12 Conclusions 

As said above the nonadiabatic calculations of dipole moments are scarce 
and few. Those of Tachikawa and Osamura [33] and Cafiero and Adamow- 
icz [35,37] are not reliable because they are not based on the molecule-fixed 
Hamiltonian operator. In our opinion the only serious attempts to the nona- 
diabatic calculation of dipole moments are those of Blinder [38,39], Kolos 
and Wolniewicz [44,52] and Wolniewicz [53,54]. All the other calculations of 
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dipole moments are based on the BO approximation. 

We have also seen that the comparison between theoretical and experi- 
mental dipole moments is not straightforward. The field-reduced line split- 
tings do not give us the square of the molecular dipole moment, as predicted 
by the oversimplified rigid-rotor model, but a kind of energy-weighted aver- 
age. Therefore, there is much to be done in this field and we hope that the 
present discussion will contribute to motivate such work. 
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Table 1: Dipole moments for some diatomic molecules 



Ref. 




Method 


LiH 


[35] 


5.8816 


non BO 


[45] 


5.879 


BO 


[33] 


6.072 


non BO 


[46] 


5.8820 (4) 


Exp. 


[26] 


5.882 ± 0.003 


Exp. 


LiD 


[35] 


5.8684 


non BO 


[33] 


6.080 


non BO 


[46] 


5.8677 (5) 


Exp. 


[26] 


5.868 ± 0.003 


Exp. 


HD 


[44] 


1.54 X 10-^ 


non BO 


[37] 


0.831 X 10-3 


non BO 


HT 


[37] 


1.111 X 10-3 


non BO 


DT 


[37] 


2.77 X 10-^ 


non BO 
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